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Abstract 

We propose a method to determine the leptonic decay constant fs* in the 
infinite heavy quark mass limit from the analysis of the radiative decay mode 
B~ — > £~&£j. The method relies on HQET symmetries and on experimental 
data from D*° -> D°-f. 
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The measurement of fs and /#» , the leptonic decay constants of the B and B* mesons, 
defined by the matrix elements: 



< 0|&7 M 75<?|£(p) >= ip„f B (1) 

< 0\h»q\B*(p, e) >= m B »f B *e^ , (2) 

represents one of the main goals of the current and future experimental investigations in 
the heavy quark physics. The reason can be found in the prime role played by fs in the 
hadronic systems containing one heavy quark. To give an example, /§ appears in the 
formula relating the mass difference between the 5°-meson mass eigenstates to \V t d\ 2 in 
the box diagram computation of the -B — B° mixing; therefore, the size of the unitarity 
triangle and the analysis of possible CP violations in the B systems crucially depend on 
the value of this hadronic parameter. 

As a second example of the relevance of the leptonic B constant, we can consider the 
Heavy Quark Effective Theory applied to the physical world of the heavy hadrons. In the 
limit rrib — ► oo, fg and scale according to the relation 

fa = is* = • (3) 

The parameter F, independent (modulo logarithmic corrections) of the heavy quark mass, 
represents a low-energy parameter related to the non-perturbative dynamics of light quark 
and gluon degrees of freedom, and plays a role analogous to in chiral theories for light 
hadrons. 

On the theoretical side, much effort has been devoted to the calculation of fa = £_ 
by non-perturbative methods such as lattice QCD ||] and QCD sum rules 0; for example, 
a QCD sum rules analysis in the infinite heavy quark mass limit provides us with the value: 

F = 0.25 - 0.45 GeVi . (4) 

(depending on the inclusion of a s corrections). 

On the experimental side, the most natural process to measure fs would be the purely 
leptonic decay mode B~ —>■ £~ve, whose decay width is given by: 



nzr r ;-,,) = ^\v ub \ 2 f B (^-) 2 ml(i-4) 2 ■ (■>) 



in Ktub ) \ m B , 



1 



If one determines \V u b\ from other processes, for example from the end-point spectrum 
of the charged lepton in the B meson inclusive semileptonic decay || |J, one can obtain 
fs by this equation. The difficulty is represented by the helicity suppression displayed 
in Eq. (|5]). Since the lepton pair must be a spin state and the antineutrino has a 
right-handed helicity, also the charged lepton is forced to be right-handed. The effect is 
the suppression factor (^-) 2 that makes the purely leptonic decay mode hardly accessible 
in the electron and in the muon case. As a matter of fact, using the value for the B~ 
lifetime r B - = 1.646 ± 0.063 ps J4], the expected rates are: 

B(B-^e-v e ) ~ 6.6 [-^1 2 [ (* ] 2 10 12 (6) 

v ; [ 0.003 J L 200 MeV J V ; 

B(B~ -> uTv u ) ~ 2.8 [-^] 2 [ tl ]2 iq-7 (7) 

v p m; [ 0.003 J L 200 Mel/ J ' V ; 

to be compared with the experimental upper bound put by CLEO 0: 

B(B~ -> e - P e ) < 1.5 10~ 5 (8) 
B[B- -> \Tv^ < 2.1 10~ 5 (9) 

(at 90 % confidence level). 

As for the channel B — > tu t , the helicity suppression is absent and the expected rate 
is larger: 

B(B- -> t~v t ) ~ 6.8 [-^] 2 [ 1 2 10~ 5 . (10) 

v ; [ 0.003 J L 200 MeV rJ V ; 

However, the r identification puts a hard experimental challenge. The upper limits found 
by CLEO § and ALEPH [§ Collaborations read: 

B{B~^t~v t ) < 2.3 10" 3 (11) 
B(B- -> t~v t ) < 1.8 10~ 3 (12) 

(at 90 % confidence level), respectively. 

For these reasons it is worthwhile to search for other paths, and analyze other possible 
decay modes that are sensitive to the value of F. For example, one could use flavour 
symmetry and consider a measurement of /#* from the spectrum of the semileptonic 
decay B — > iriu^ near zero recoil 0, || compared to the spectrum of D — > 7r£z/^ in the same 
kinematics. In this case, however, one has to face a strong phase-space suppression. 



Another possibility, suggested in refs.Q and flOR , is provided by the radiative leptonic 
decay channel 

B~ -> H~v^i (13) 
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which does not suffer of the helicity suppression because of the presence of the photon 
in the final state. In this case there are several uncertainties related to the hadronic 
parameters appearing in the matrix element governing the decay mode (0). Within 
such uncertainties, the branching ratio of the decay fljJD has been estimated in the range 
10~ 7 — 10~ 6 (in the case of light leptons), thus making the channel promising for a future 
£>-factory. We shall now study how this decay mode can be used to measure fs* ■ 

In order to analyze the dependence of the amplitude for the process (|T3"D on f B *, let 
us consider the diagrams which describe it; they can be divided into two classes. The first 
class consists of structure dependent (SD) diagrams such as those of Figs. 1; the second 
class contains bremsstrahalung diagrams where the photon is emitted from the B~ or 
from the charged lepton leg. 

The bremsstrahalung amplitude is given by: 

Mb = l f B e^V ub m,[(F(k^ - ZJ*)ji(l - l5 )v - ^—^4 $ (1 - 75 H (14) 

where p, pi, k are the momenta of B~ , /z and 7, respectively , e is the photon polarization 
vector, and F(k 2 ) is the B~ electromagnetic form factor. This contribution vanishes in 
the limit —>■ and we shall make this approximation, so that the relevant diagrams 



governing (p~3|) are the SD diagrams. We shall suppose, as in fllOfl , that in these polar 
diagrams the intermediate state can be a J p = 1~(B*) and a positive parity B** meson. 
The amplitude with intermediate P(= B*,B**) state can be written as follows: 

MU = ^V ub A(B - Pi) j- _ JL _ ^ < 0|«7*(1 - !,)b\P > I, , (15) 

where = £(p/)7 At (l — 7s)^(^) is the lepton current, A(B — > P7) is the amplitude of 
the process B — > P7, and P indicates the pole. From Eq.(|I^) (with P = B*) it is clear 
that, for light leptons in the final state, the radiative B decay can give access to the decay 
constant /#* provided that 

1) one has enough information on the amplitude A(B* — > -B7), 

2) the remaining part of the SD contribution J\z[g£) is small compared to J^Agjj . 

Let us begin by discussing the first point and let us consider the contribution of the 
B* pole to (|15D . In computing the on-shell amplitude A{B* — > B^) one has to take into 
account the coupling of the electromagnetic current to the heavy quark and to the light 
quark, i.e. the terms arising from the decomposition J^ m = 6567^6+ e g g7 M g (ef>, e q = quark 
charges) . 



The b— quark contribution is described (in the limit mj — > oo) by the amplitude: 

e*» < B{v')\e h b ltl b\B*{v : r,) >= i e b £(y ■ ^ V u v a v ,/3 , (16) 

where v and v' are B* and B four-velocities, respectively, and £(v ■ v') is the universal 
Isgur-Wise form factor (IW) For on-shell B and B*, since v ■ v ' = ^^tragl — ^> we 
can use the normalization of the IW function: £(1) = 1. 

The second contribution: e*^ < B\e q q^^q\B* > represents the coupling of the electro- 
magnetic current to the light quark q (q = u for B~ decay); this contribution dominates 
in the m;, — > oo limit, and is more uncertain since it cannot be estimated within HQET 
because it involves light degrees of freedom. On the experimental side we have no data 
on the width B* — > B*y at the moment, and it is unlikely it will be measured in the near 
future. On the other hand, the experimental D* 0, ~ — > D°~ , -f branching ratios are known 
(even though the full D*°~ width has not been measured yet) and we may presume that 
in future we will get information on the D* partial radiative width. 

Our proposal is to use these pieces of information to obtain B* — > Bj. 

The D*° radiative width is given by: 

r(/r° -> d^) = A^B±^ (17) 

127T mpo ' 



where g 7 is the photon momentum in the D* rest frame and 



9d*d-j — e 



m D * A 



q 



(18) 



(e c = 2/3, e q = e u = 2/3). A measurement of T(D*° — > D°j) would provide a determina- 
tion of the mass constant A q that parametrizes the matrix element < D^'y^D* >. On 
the other hand 



A(B*(v,r))^B(v')y(q,e)) 



eb e q 



x 



_m B , A 9J 
m B *y/m B mB*e llva pe*^rfv 0l v' p , (19) 



and therefore, by the the knowledge of A q , one would obtain the matrix element needed to 
compute Eq.flloT) (with P = B*). In other words with the approach we have just described, 
we might extract F by the measurement of both T(D*° — > D*j) and &(B~ — > ^P^j); we 
shall discuss below the sensitivity of this approach. 
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In order to assess the reliability of this method we shall also discuss as a consistency 
test based on information in part already available, i.e. data onS-> tt£u [12], D — > Klv 



13| and on the BR's B(D* -> Dn) and B(D* -> D7) [0. Let us consider the partial 



width 

r(zr° - dV) = ( 9 -Ep\ 2 (20 ) 

where q n is the pion momentum in the D* rest frame and Qwd-k is the strong D*Dir 
coupling constant, which, in the heavy quark mass limit, is given by 

9D*Dn = 9d*+d^+ = ^j-y/rn D m D , . (21) 



The ratio of the partial widths (|20|) and ( fTF]) is experimentally known [13] 



from which the ratio of the Z)* decay constants can be obtained: 

= 0.808 VR GeV . (23) 

In order to extract go*D-y from Eq.(p3|) one needs the strong coupling constant g; several 



determinations of g have appeared in the literature, based on QCD sum rules M, 14, 15 



or potential models 16|| ; they indicate a value in the range 0.2 — 0.4. We shall avoid, 



however, to rely on these estimates and we shall try to minimize the theoretical bias using 
experimental data of the decay D -> Kiv (B(D + -> K°e + u e ) = 6.6 ± 0.9 10~ 2 [|TJ) and 
the recently observed decays [[LJ 



B(B^ir£v) = 1.70 ± 0.50 10" 4 (BSW) (24) 

B{B^Ttiv) = 1.19 ± 0.65 10" 4 (ISGW) . (25) 

The two determinations refer to the model used in the Montecarlo code to compute the 



efficiencies: the Bauer et al. model [17] or the Isgur et al. model [18 



For these decays, assuming a simple pole formula for the form factors f+~* n (q 2 ) and 
f+^ K {q 2 ), which is what is experimentally found for D + — > K°£ + U£ and generally accepted 
on theoretical grounds [|19|, ^(J, |21], |22|] for B — > ntv, one gets the following results. The 
semileptonic B — ► -niv partial width is given by 

r^^HFVi^^/W (26) 



with 



J 



Q 2 

I m ax 



dq 2 



m B ,ml,q 2 ) 



A 3 / 2 (- 2 
(l-g 2 /m|*) 2 



In ( P7|) A is the triangular function. For D — > If^i/g one has to change mg 
Kfe -> Kg, m B * -> mi?., — > mjf. From Eq.(g7J) one gets: 



(27) 



2.95 10" 



II' 



«6 



1/2 



Ge^i = (1.07 - 1.28) 10" 1 



1 



GeV 3/2 , (28) 



I 0.003 I 



Fg 



0.45 



1/2 



Geyi = 1.18 VQ~ x GeV*l 2 \ 



(29) 



thus, the use of B — > 7r£z/£ data or D — > Kev e , employing light and heavy flavour symme- 
tries, give similar results. Eqs. (|28|j29| ) allow to obtain, from (|3|), Qd*d^ and therefore, 



from (|i~8D, A g as a function of F. 



We shall use below this approach as a consistency test of our method. In both cases, 
the proposed method and the consistency test, after having put the value of A q into ([03), 
we can compute A(B* — > S7) and the amplitude , i.e. the contribution of the 



B* pole to the decay B 
follows: 



( B* ) 

The final expression for the amplitude Ai SD is as 



M 



(B*) 



Cxf. 



B* 



SD (vk + A) 
where A = mg. — tub, and C\ is given by: 



Ci = —H= V ub J-^- y/m B m B * e 
V2 2m B 



m B * 



From ( PH| ) we can compute the contribution of the B* pole to B[B 
function of the parameter F. 



(30) 
(31) 

fjTi/fi'y) as a 



Let us study the effect of the B** pole. It is well known that in the limit of infinitely 
massive heavy quarks (tuq — > 00) the strong dynamics of the heavy quark decouples from 
the dynamics of light degrees of freedom, with the consequence that the spin of the heavy 
quark and the spin of light degrees of freedom are separately conserved. In the case of 
the first orbitally excited heavy states, with orbital angular momentum I = 1, the total 
angular momentum of the light degrees of freedom can be sg = 1/2 or = 3/2. To 
each of these values corresponds a doublet of positive parity heavy mesons: the doublet 
(B , B[) with J p = (0 + , 1 + ) in corrispondence to si = 1/2, and the doublet (B 1} B 2 ) with 
J p = (l+,2+) (s e = 3/2). 
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In the charm sector only the members of the se = 3/2 doublet (both in the strange 
and non-strange charmed channels) have been observed: in the non-strange case they are 
the states _Di(2420) and .02(2460) decaying to Dtt, D*tt with pions in .D-wave, therefore 
with small decay widths (T Dl = 18^ MeV and V D2 = 23 ± 10 MeV flip . The states 



se = 1/2 can decay into Dtt, Dtt* with pions in 5*-wave, and their large decay width 
makes their observation rather problematic. In the beauty sector, recent results from 
LEP Collaborations show the existence of positive parity orbitally excited B mesons 
(B**) with an average mass m B » = 5732 ± 5 ± 20 MeV |23| . 



Only the axial vector states 1 + can contribute as poles to the decay (Fig. 1 b); 
moreover, the state Bi, having si = 3/2, has vanishing coupling to the weak current in 



the limit m& — > oo [p2], and only the state B[ (with si = 1/2) gives a contribution in the 
same limit. 



As to B[, following the same steps leading to Eq.(|30|) we get: 

M ™ = ' (v^k+'W 6 " v *" ~ v ' * e ^ )P (32) 

where A' = tub' — ttib , 

< ^\h^Q.\B[(P^) >= fBtjniXjIn , (33) 

and 



G F m B > r 2e b r 1/2 (l) e q 



(34) 



The function Ti/ 2 {v ■ v') is the universal form factor, analogous to the IW function, de- 
scribing the matrix element of the weak current between positive parity heavy mesons 
and the doublet (B, B*) in the limit mj, — > oo [E3]; it can be defined as follows: 



< B[(v', r])\b-f^ 5 b\B(v) >= i t 1/2 (v ■ v')^mBmB' x e^p rf v v a v' p , (35) 

and its value at v ■ v' — 1 has been estimated by QCD sum rules ||25|| : T!/ 2 (l) — 0.24. We 
notice that, because of the factor l/mg and the small value of 71/2(1), the first term in 
the r.h.s. of (|35| ) is expected to be small as compared to the second one. 



Putting together the contributions of B* and B[, we get: 

12 f2 



2 rm B /2 



{E 1 + A) 2 (E 7 + A') 

r (B«) + V {B[) (36) 



The relative contribution of the B* and of the B[ poles is given by: 



r(Bj) = psw \c 2 f B[ \ 2 

T(B*) PS( B *) ICJb*] 2 { ' 

therefore it is weighted by the ratio of the phase-space coefficients 

PS( B D 

j^m * °- 54 ^ 

(using the experimental values of the mass differences A = 46 MeV and A' ~ 500 MeV) 
and depends on the ratio (C2f B > 1 )/(Cif B *). This last ratio is basically determined by the 
coupling of the electromagnetic current to the light quarks. This can be seen from Eqs. 
fl3i"| , |34]) neglecting l/m b terms, since the heavy quark magnetic momentum is subleading 
in the inverse heavy quark mass expansion: 



C 2 f B { „, m B{a fs[ A g 

Ci f B * m B * f B * A' 



(39) 



The ratio can be estimated using experimental results from semileptonic D— meson 



decay channels |22| |26fl. As a matter of fact, one can make use of the Vector Meson Dom- 
inance (VMD) |26| ansatz and assume that the coupling of the electromagnetic current 
to the light quarks is mediated by vector meson states V (V — p, u). By this assumption 
the matrix elements < 5|g7 M g|S* > and < B\q , ~f >l q\B' 1 > can be written as: 



< B(p')\q^q\P(p,e 2 ) >= £ < B{p')V{q, e 1 (rj))\P(p, e 2 ) > 



i 



v, v f - m v 

<0\Wtf\V(q,e 1 (rj))>, (40) 

where P = B* , B[ and therefore they are proportional to the strong couplings < B V\B* > 
and < B V\B[ >; these matrix elements have been estimated in [^2| in the framework 



of an effective chiral theory for heavy mesons, using experimental information from the 
D — > K*lv semileptonic decay and the (approximate) symmetries of the effective theory 
(heavy flavour symmetry and flavour SU(3) symmetry for the light degrees of freedom): 

< B{v')V{k,e)\B*{v,r}) >= e^ v v a v'^m B m B .m B ^g y \ (41) 

< B(v')V(k, e)\B[(v, rj) >= [(k ■ v)(e* ■ rj) - (k ■ rj)(e* ■ v)]^m B m B ,2V2g v p (42) 
where A = -0.4 GeV' 1 , p = -0.13 GeV~ x and g v = 5.8 J22jj . As a result one finds 

fw A 

^^0.4 . (43) 

JB* Ag 
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Therefore, we conclude that 



iw - o- 1 ( 44 ) 

i.e. the B[ pole represents a contribution at the level of 10% of the radiative B~ — > ^v^i 
width. The phase-space suppression is even larger for the contribution of other (radially or 
orbitally) excited B states, which therefore can be safely neglected. Then, we can assume 
that the diagram of Fig. lb adds a term of the order of 10% as compared to Fig. la; 
therefore we shall take into account it by multiplying the contribution of the B* pole to the 
width by a normalization factor K — 1.1. We are aware that this estimate of the diagrams 
of Fig. lb is more uncertain than the evaluation of Fig. la. Nevertheless we observe that 
the VMD hypothesis successfully describes a number of low energy phenomena involving 
photon radiation |27[] ; moreover these theoretical uncertainties will not affect strongly the 
final determination of B{B — > fii"~f), since, in any case, the contribution of Fig. lb is 
much smaller than the B* term. 



The main difference of the above analysis with respect to ref. (10| consists in the 
evaluation of the coupling g BB j 7 and in the role of Qb*b^- In [0 non relativistic quark 
model has been employed to compute gB*B-y, using A q = 330 MeV ~ m u (constituent 
mass of the light quarks); moreover, it has also been adopted to estimate Qb' x b^ — 3i ^| 7 ■ 
On the other hand, in our approach gs'B-y is a quantity that should be inferred from 
experimental data, thus reducing the dependence of the results on the hadronization 
model. The prediction in [10] for BiB — > /iz/7) is in the range 10~ 7 — 10~ 6 . 



Let us now come to the numerical results. As we observed above, to set definite predic- 
tions one needs the experimental input T(D*° — > D 0/ y), which is not available yet. In order 
to test the sensitivity of the method, we have used two theoretical determinations of the 
radiative D* width: T(D*° -»■ L>° 7 ) = 22 KeV (ref.0) and T(D*° -»■ £>° 7 ) = 11 KeV 
(ref. p8|| ) ; with these input data, which represent rather extreme cases (an intermediate 



prediction can be found in [p6fl ), and using V u b = 0.003, we get the curves (a) and (b) 
depicted in Fig. 2. From this figure one can see that the branching ratio of the radiative 
decay mode B~ — > \i~v^j is expected to be larger than in the purely leptonic mode. This 
can be appreciated by comparing the prediction in Eq.(7) with the outcome in Fig. 2 in 
correspondence to the value F ~ 0.45 GeV 3 ^ 2 : B(B~ — > yTv^) = 4 — 10 x 10~ 7 (de- 
pending on the D*° D 0, y decay width). This implies that the radiative mode, for quite 
large values of the leptonic constant F is favoured by a factor 2 — 3 (at least) with respect 
to the purely leptonic mode B~ — > \Tv\ moreover, as far as the statistics is concerned, 
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the use of radiative decay allows to employ the electron channel as well, which provides 
a gain of an additional factor of two. 



To test this result using the semileptonic data (Eqs. j29|) , we have reported in Fig.3 
two curves corresponding to the central values for B — > nlv given in Eq. ([24]) and (p5j); 
for the range of F in Eq.([|) we obtain &{B~ — > ^Tv^f) ~ 2 — 3 x 10~ 7 . The data for 
D — > Klv give an intermediate result. As we have stressed above, the outcome in Fig.3 
is based not only on experimental data, but also on additional theoretical assumptions, 
such as the polar dependence of the semileptonic B — > tt form factor. So, such results 
are characterized by a quite large uncertainty; nevertheless, for intermediate values of F 
there is an overlap region with the result in Fig. 2 that allows us to conclude that the 
consistency test does not contradict our main results in Fig. 2 

In conclusion, our analysis confirms that the decay channel B~ — > y^v^ can be used 
as a way to measure the leptonic B* decay constant; one expects larger decay rates than 
in the purely leptonic case, even though the detection of the photon in the final state 
may reduce the reconstruction efficiency Moreover, since the method described here 
is strongly based on HQET symmetries, it would be useful in any case to look for the 
radiative leptonic B decay to test experimentally HQET predictions in this context. 



Note added 



After completing this work we became aware of the paper |29| where the calculation of 
B — > ivy is performed by light-cone sum rules. The result confirms the expectation that 
the radiative decay rates are larger than the purely leptonic rates. 
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Figure Captions 



Figure 1 

Diagrams dominating the B~ — > t'vf) decay mode in the limit mg — > 0. B* is the vector 
meson state, B** is the 1 + axial vector meson state. 

Figure 2 

Branching ratio B(B~ — > ^v^) obtained according to the method described in the text. 
The curves (a) and (6) refer to the values: T(D*° — > D°7) = 22 ifeV |TBj (continuous 
line) and r(D*° -> D° 7 ) = 11 ifeV || (dashed line). 

Figure 3 

Branching ratio B(B~ — > //"tyy) obtained using data on the semileptonic 5 — > 7r£z/ decay. 
The curves (a) and (6) refer to the input values: £>(!? — > niu) = 1.70 10~ 4 (continuous 
line) and B(B -> vr^) = 1.19 1(T 4 (dashed line). 
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Fig. 1 a) 




Fig. 1 b) 
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fig.2 



